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Abstract. Over the last few years, the existence of mutual feedback effects between accreting 
supermassive black holes powering AGN and star formation in their host galaxies has become 
evident. This means that the formation and the evolution of AGN and galaxies should be con- 
sidered as one and the same problem. As a consequence, the search for, and the characterization 
of the evolutive and physical properties of AGN over a large redshift interval is a key topic 
of present research in the field of observational cosmology. Significant advances have been ob- 
tained in the last ten years thanks to the sizable number of XMM-Newton and Chandra surveys, 
complemented by multiwavelength follow-up programs. I will present some of the recent results 
and the ongoing efforts (mostly from the COSMOS and CDFS surveys) aimed at obtaining a 
complete census of accreting Black Holes in the Universe, and a characterization of the host 
galaxies properties. 
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Since the launch in 1999 of both the XMM-Newton and Chandra satellites, a large 
(> 30) number of surveys covering a wide fraction of the area vs. depth plane (see 
Fig. 1 in |Brandt fc Hasi ngcr 2005] see also IHickox 2009[) have been performed and our 
understanding of AGN properties and evolution has received a major boost. Thanks 
to vigorous programs of multiwavelength follow-up campaigns which, since a few years, 
have became customary, sensitive X-ray observations turned out to be highly efficient in 
unveiling weak and/or elusive accreting black holes, in a variety of otherwise non- active 
galaxies, such as (among others): X-ray Bright Optically Normal Galaxies (XBONG, 
IComastri et al. 2002[) . Extremely Red Objects (e.g. IBrusa et al. 2005|) . Sub Millimeter 
Galaxies (e.g. Alexander et al. 2005) high-z starforming systems (e.g. ?. IFiore eFa l. 2008). 
In many of these cases, the AGN responsible for the X-ray emission is overwhelmed at 
longer wavelength by the host galaxy light and/or the obscuration might be connected 
to processes within the host galaxy itself, such as the star formation rate and the pres- 
ence of dust lanes or starburst disks (see e.g. |Hopkins et al. 2009[ ). This suggests that the 
accretion activity (especially in high-redshift sources) is unambiguously revealed thanks 
to the presence of a strong X-ray emission (see e.g. discussion in IBrusa et al. 2009a[) and 
therefore, the combination of both X-ray and optical classifications can be crucial to fully 
assess the nature of the candidate AGN. 

The high level of completeness in redshift determination for a large number of X-ray 
selected AGN (up to a few thousands) has made possible a robust determination of the 
luminosity function and evolution of unobscured and mildly obscured AGN which turned 
out to be luminosity dependent: the space density of bright QSOs (Lx > 10 44 erg s _1 ) 
peaks at z~ 2-3, to be compared with the z~0.7-l peak of lower luminosity Seyfert galax- 
ies (|Ueda et al. 20031 lETFranca et al. 20051 ISilverman et al. 20"08al lEbrero et al. 2009[) . 
Based on these works, Marconi et al. (2004) and Merloni (2004) were the first to pro- 
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pose that SMBH undergo a "anti-hierarchical" evolution, in the form a differential 
growth (earlier and faster for larger black holes) . This anti-hierarchical behavior observed 
in AGN evolution (similar to that observed in normal galaxies, e.g. ICowie et al. 1996|) 
provided an important, independent, confirmation that the formation and evolution of 
SMBH and their host galaxies are likely different aspects of the same astrophysical prob- 
lem. In particular, the accretion onto supermassive Black Holes triggered by galaxies 
mergers and/or collisions may provide the "feedback" needed in almost all recent models 
for galaxies formation and evolution, in order to recover the properties (masses, lumi- 
noisity, clusetring) we observe today in local systems. 

In this general framework the differences between "obscured" and "unobscured" AGN 
are no longer and uniquely described under a zeroth-order geometrical unification model 
(in which they are simply related to orientation effects. lAntonucci fc Miller 1985)) . but 
can be interpreted as due to the fact that the same object is observed in different evolutive 
phases. This hypothesis is further supported by the finding that absorption is much more 
common at low luminosities (jUeda et al. 2 003, Ma iolino et al. 2007[1 and at high redshift 
(|La Franca et al.~20 05. Treister & Urry 2006, Hasinger 2008) as emerges mainly from X- 
ray surveys. The luminosity and redshift dependence of the obscuring fraction may be 
naturally linked to the AGN radiative power (related to the intrinsic X-ray luminosity) 
which is able to ionize and expel gas and dust (more common at high-z) from the nuclear 
regions, nicely fitting the current framework of AGN formation and evolution sketched 
above (see, for examle, [Hopkins et al. 2006| and this volume). 

The complete picture is likely to be much more complicated, depending on many other 
parameters (such as, e.g., the BH mass, the Eddington ratio, the QSO duty cycle), and in 
particular related to the complex "lightcurves" of AGN, which, in turns, depend on the 
detailed hydrodynamics adopted in the simulations. A correct and complete identification 
of unobscured, obscured and highly obscured AGN at all redshifts (and especially in the 
z=l-3 interval, where most of the feedback is expected to happen) is therefore crucial for 
a comprehensive understanding of the still little explored phase of the common growth 
of SMBHs and their host galaxies. 

Here I highlight recent results the characterization of the high luminosity tail of the 
X-ray luminosity function from the COSMOS survey, and on the host galaxy properties 
of obscured AGN at cosmological distance, from ananalysis of X-ray selected sources in 
the Chandra Deep Field South (CDFS). I also present evidence for an evolution of the 
Mflff-Mfluije relation as derived from optically selected BL AGN from the zCOSMOS 
survey (almost all of them also X-ray emitters). 

1. The high-luminosity tail of X-ray selected AGN: results from the 
COSMOS survey 

The high completeness in optical identifications of the XMM-COSMOS sample ( Brus a et al. 20 09b) 
and the availability of spectroscopic and photometric redshifts ( Sal vat o et al. 2009]) al- 
lows us to estimate the number densities of AGNs as a function of luminosity and redshift. 
In calculating the binned number density evolution we considered only the subsample of 
sources detected in the hard (2-10 keV) band, an included in the flux limited sample (925 
sources at S x = 6 X 10~ 15 erg s _1 cm~ 2 ; ~ 60% with spectroscopic redshifts available). 
Since the hardness ratios between the hard and soft bands give the estimates of the ab- 
sorption column density, we present the XLF of the de-absorbed 2-10 keV luminosity. In 
our first analysis, we use the traditional ^V^a estimator (Avni&Bahcall 1980). Addi- 
tional details on the method used in deriving the XLF are given in Brusa et al. (2010) 
and Miyaji et al. (in preparation). 
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We compared our first estimates of the number density evolution of the most lumi- 
nous AGN population with recent measurements. The results are presented in Figure [TJ 
where the number density of XMM-COSMOS sources is plotted versus the redshift in 
the luminosity bin LogLj^ = 44.5 — 47 (squares). We only included the bins where the 
lowest L x lowest z boundary corresponds to a 2-10 keV flux above our limit for plot- 
ting. In the same figure we plot the data points obtained from the AEGIS survey by 
Aird et al. (2010)] circles), along with their best fit luminosity function (black curve). 
The results based on 6 datapoints within z=l-3, with ~ 20 objects each may provide ro- 
bust constraints on the shape of the XLF. In particular, the XMM-COSMOS data points 
seem to favour a higher redshift (z~ 2) peak for the space density of luminous quasars 
as predicted by a LDDE parameterization, represented in the figure by the Ebrero et 
al. (2009, short-dashed curve) and the Silverman et al. (2008, long-dashed curve) lumi- 
nosity functions, though different in the details, and reproduced by XRB synthesis mod- 
els ( |Gilli, Comastri fc Ha singer 2007; dot-dashed curve) rather than the lower redshift 
(z~ 1) peak expected from LADE model recently proposed by Aird et al. (2010). 

2. Black Hole and Star formation activity at z> 1 

In order to study the host-galaxies of obscured AGN, we defined a sample of 116 "bona 
fide" obscured AGN, detected in the lMs CDFS observation (jAlexander et al. 2003|) and 
for which deep Infrared (IRAC and MIPS), K-band and multiwavelength optical pho- 
tometry is available (from IGrazian et al. 20061 the MUSIC survey) We selected sources 
without broad lines in the optical spectra and with small optical nuclear emission with 
respect to the host galaxy optical emission. 
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Figure 1. The number density evolution from the XMM-COSMOS 53 field AGN (squares) 
compared with recent results from Aird et al. (2010; circles, and black solid line), Ebrero et al. 
(2009; short-dashed curve), Silverman et al. (2008; long-dashed curve), and the expectation of 
XRB synthesis models taken from Gilli et al. (2007, dot-dashed curve). 
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We found that the hosts of obscured AGN are redder in (U-V) rest frame than the 
overall galaxy population at the same redshift: in particular, obscured AGN mainly 
populate the red sequence and the green valley in the color-magnitude plots in agreement 



with the results of Silverman et al. (2008b) For the MUSIC sample the U-V galaxy colors 
are strongly correlated with the K band absolute magnitude (Figure [5]), and therefore 
with the galaxy stellar mass, with the most massive systems having a redder color. The 
hosts of the obscured AGN are therefore found in the red-massive tail of the distribution 
of optically selected galaxies in all three redshift bins considered (see labels in the figure). 
AGN feedback is often invoked as one of the main responsible for the observed galaxy 
colors (e.g. Nandra et al. 2007, Hasinger 2008). However, it is well known that the main 
ingredient for nuclear activity is the presence of a SMBH in galaxy nuclei, and that 
SMBHs are found nearly exclusively in massive galaxies (e.g. Magorrian et al. 1998). 



Therefore, it is not truly surprising to find AGN hosted in massive galaxies and the 
simple presence of AGN in massive red galaxies is not enough to argue for a significant 
feedback effect on the observed colors, because of the strong color-mass correlation. Were 
AGN feedback responsible for the observed red colors, since galaxy colors are strongly 
correlated with the galaxy mass and AGN are found preferably in massive galaxies, then 
AGN feedback should be also considered as one of the main players in the building of 
the galaxy mass-color correlation. 

We found that about 2 /3 of the obscured AGN hosts at all redshifts show substantial 
(> 10 M© yr -1 ) star formation activity (Figure El where two histograms of the SFR 
derivd from SED fitting are reported, for optical and Infrared selected smaples) and 
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Figure 2. U-V(rest frame) vs. the absolute K-band magnitude for the obscured AGN (squares 
and triangles) and the underlying galaxy population (contours). The upper panels show the 
comparison of the X-ray selected sources with the optically-selected (z-band) galaxy population 
in three different redshift bins: z=0. 6-1.0, z=l-2, and z—2A from left to right, respectively. The 
lower panels show the comparison for the subsample of objects detected also at 24 micron, in 
the same redshift bins. 
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about half live in galaxies which are still actively forming stars with respect to their 
mass. For these sources, the observed red colors are likely due to dust extinction rather 
than evolved stellar population. We then conclude that a significant fraction of obscured 
AGN live in massive, dusty star-forming galaxies with red optical colors. 

We compared the number of obscured AGN and of all X-ray selected AGN to the 
number of field galaxies in broad bins of galaxy stellar mass (M* = 10 10 — 10 12 M Q ) and 
redshifts (z=0.6-l, z=l-2, z=2-4). We find that the AGN fraction increases with the host 
galaxy stellar mass, from ~1% at M» ~ 10 10 M to -30% at M» ~ 3 x 10 11 M (see 
also lYamada et al. 2009), and the actual trend of increasing AGN fraction as a function 
of the stellar mass is probably steeper given the uncompleteness of the MUSIC sample at 
M* < 10 11 M© (see Figure HJ. We compared this trend with that observed in the local 
Universe ([Best et al. 2005[) for AGN with luminosity above similar thresholds. While the 
observed trend is the same, in all the investigated redshift bins the AGN fraction is higher 
than that observed in the local Universe, and it could likely be even higher. In fact, we 
are comparing here AGN selected with two different methods: forbidden emission line lu- 
minosity (SDSS) and X-ray emission (GOODS). The latter sample does not contain most 
Compton thick AGN. On the other hand, Compton thick AGN may well be present in 
[OIII] selected AGN samples. The fraction of Compton thick AGN not directly detected 
in deep Chandra surveys is estimated between 40% and 100% of the X-ray selected AGN, 
using infrared selection or other techniques (see Donley et al. 2008 IFiore et al. 20081 and 
references therein). Therefore, under the simplest assumption that this Compton Thick 
AGN fraction is constant with the galaxy mass, the discrepancy observed in Figure 10 
can increase by up of a factor of 2. 

The fraction of active galaxies to the total galaxy population is proportional to the AGN 
duty cycle. Our results would thus suggest higher AGN duty cycles at z=2-4 than at 
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Figure 3. The distribution of the SFR for both the optically (upper) and MIPS (lower) selected 
sample of field galaxies (open histogram) and obscured AGN (filled histogram) in the redshift 
interval z=l-4 and in a specific stellar mass range (M* = 10 10 — 10 12 Mq). 
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z=0, in agreement with expectations from most recent semi-analytic models (e.g. ?), in 
which at higher redshift the AGN activity is present in a larger number of galaxies than 
locally. 

The fact that the most luminous obscured AGN are found in the most massive galaxies 
at all investigated redshifts may suggest that the L/L^dd of the obscured AGN is similar, 
particularly in the case of the most luminous sources (logLx> 43 erg s _1 ), for which the 
threshold in luminosity introduces a bias against the sources accreting at lower rates in 
the lowest redshift bin. Assuming the local Magorrian relation between Mbh and M* (e.g 
Marconi & Hunt 2003) and a bolometric correction of 20 (e.g. Marconi et al. 2004) the 
median observed values of Lx/M* correspond to L/LEdd ~ 0.1. Although suffering from 
large uncertainties associated with the stellar mass and BH mass estimates, this value 
can be considered as representative of the accretion state of the most luminous, obscured 
AGN in the present sample. Similar results are obtained for Chandra Deep Field North 
X-ray sources at z=2-4 (Yamada et al. 2009P and are also typical of unobscured Type 1 
AGN at z> 1 (see Merloni et al. 2009). 



3. The evolution of ~M.BH-^-Buige relation 

Local scaling relations have proved themselves unable to unambiguously determine the 
physical nature of the SMBH-galaxy coupling. A large number of theoretical models for 
the AGN-galaxy interaction have been proposed, all tuned to reproduce the z=0 observa- 
tions. One obvious way out of this impasse is the study of the scaling relations' evolution, 
which had, until recently, been limited to a handful of objects in narrow redshift windows 
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Figure 4. Fraction of obscured AGN with Lx > 10 43 erg s _1 as a function of the stellar mass 
in three different redshift bins. Filled circles refer to the "obscured" AGN sample. The dashed 
(continuous) lines represent the fraction of AGN with L(OIII)> 10 (10 ) Lq in an optically 
selected sample in the local Universe (z< 0.2) from the SDSS (Best et al. 2005}. 
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(Peng et al. 2006j & this volume; IWoo et al. 20081 & this volume, IJahnke et al. 20091 & 



this volume). 

Within the zCOSMOS collaboration (a massive spectroscopic campaigns with ~ 10.000 
spectra of galaxies and AGN, [Lilly et al. 2007[ ), we have pioneered a new method to unveil 
the intrinsic physical properties of AGN hosts, using the unprecedented multi- wavelength 
coverage of the COSMOS field. In |Merloni et al. (2009)| we have been able to measure 



rest frame K-band luminosities and total stellar masses, A/», of the hosts of 89 broad-line 
(un-obscured) AGN in the redshift range 1 < z < 2.2, for which we measured the black 
hole mass, Mbh, through the virial method (see e.g. lPeterson et al. 20 04 & this volume). 
This sample constitutes the largest high redshift sample so far for which reliable black 
hole and galaxy masses are available. We found that, as compared to the local value, the 
average black hole to host galaxy mass ratio appears to evolve positively with redshift, 
with a best fit evolution of the form (1 + 2 )°-68±o.i2 ^ gee pjg urc [5])_ A thorough analysis 
of observational biases induced by intrinsic scatter in the scaling relations reinforces the 
conclusion that an evolution of the Mbh — Af* relation must ensue for actively growing 
black holes at early times: either its overall normalization, or its intrinsic scatter (or 
both) appear to increase with redshift. 
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Figure 5. Redshift evolution of the offset measured for zCOSMOS type 1 AGN form the local 
Mfljj-M, relation. Different symbols identify different ranges of Eddington ratios with upwards 
arrows representing upper limit on the host galaxy mass. The solid black line shows the best 
fit obtained assuming an evolution of the form ALog(MBH /Af*)(z)=<J(l+z) with <5=0.68±0.12. 
The solid, dashed and dot-dashed lines show the bias due to an intrinisc scatter in the scaling 
relation of 0.3, 0.5 and 0.7dex, respectively. In the inset we show a comparison of our data (black 
circles, binned) with data collected from the literature. 
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